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We study the derivative expansion for the effective action in the framework of the Exact Renor- 
malization Group for a single component scalar theory. By truncating the expansion to the first two 
terms, the potential Uk and the kinetic coefficient Zk, our analysis suggests that a set of coupled 
differential equations for these two functions can be established under certain smoothness conditions 
for the background field and that sharp and smooth cut-off give the same result. In addition we 
find that, differently from the case of the potential, a further expansion is needed to obtain the 
differential equation for Zk, according to the relative weight between the kinetic and the potential 
terms. As a result, two different approximations to the Zk equation are obtained. Finally a numer- 
ical analysis of the coupled equations for Uk and Zk is performed at the non-gaussian fixed point in 
D < 4 dimensions to determine the anomalous dimension of the field. 
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I. INTRODUCTION 



The dynamics of a quantum system with very many (an infinite number of) degrees of freedom, the subject of 
Quantum Field Theory, involves very many scales at once. As it is by now well understood this is the source 
of the main difficulties encountered in the study of such a system. In perturbation theory one of the most striking 
manifestation of these difficulties is the appearance of divergent quantities. The origin of this nonsensical result should 
probably be addressed to a wrong choice of the representation. In fact, unlike any problem with a finite number of 
degrees of freedom, in the present case the choice of the representation is part of the solution of the dynamical problem 
due to the existence of inequivalent representations, but up to now not very much progress has been made in this 
direction. Moreover there are phenomena, such as the chiral symmetry breaking, where it is necessary to resort to 
non perturbative methods, since no small parameter is available for a perturbative expansion. Typical tools for this 
kind of investigations are the Schwinger-Dyson equations as well as the variational approximations. Unfortunately 
even within these approaches divergent quantities are still present and, unlike the perturbative case, a systematic way 
of handling these divergences is still lacking. This has cast serious doubts on the validity of the results obtained with 
these methods. 

From a different stand point the Renormalization Group Method, more precisely the Wilson formulation of the 
Renormalization Group (RG) M, has been invented to handle those dynamical problems involving very many scales. 
The main idea of the method is to break the original problem into a series of subproblems each involving a more 
restricted scale range and to seek for its solution by iteratively solving these subproblems. The advantage of this 
method is that it is non perturbative by construction and that no problem of divergences is ever encountered. 

There are different (more or less equivalent) ways to implement this program. The realization in terms of differential 
equations in the momentum space generally named 'exact Renormalization Group equations', has been widely 

studied in the last years and even employed as a tool for investigating many problems related to phase transitions 
where the flow of the physical parameters cannot be controlled by simple perturbative calculations. 

Among the various existing formulations of the exact RG equations, the Wegner-Houghton (WH) equation j| is 
particularly appealing for the transparency of its physical interpretation, and the most straightforward approximation 
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to this equation, obtained by restricting the problem to analyzing the RG evolution of the costant (x-independent) 
part of the action, namely the Local Potential Approximation (LPA), has been shown to contain many interesting 
features. For instance the exactness of the LPA for N — > oo has been found in Q (N is the number of scalar fields 
<f>i), as well as the fixed point structure of the scalar theory has been recovered in the LPA furthermore, the 

typically non-perturbative feature of convexity of the effective potential is already contained in the LPA ■ 

In order to extract informations on the dynamics of the field, however, it is necessary to improve the LPA by 
including the effects of the field fluctuations. The natural extension of the LPA is the obtained by considering a 
Derivative Expansion of the action, that is by adding in the general form of the action terms with higher and higher 
number of derivatives of the fields. This expansion is motivated by the expectation that, for a sufficiently smooth 
background configuration, the action should be quasi-local. A set of two coupled equations for the local potential 
Uk(4>) (k is the momentum scale at which the action Sk is defined) and the coefficient Zk{4>) of the lowest order 
derivative term, d^d^tfi has been deduced and studied by resorting to a smooth cut-off [||,|llj that allows a weighted 
integration of the various modes appearing in the original action. An equation for Zk(4>) which allows to recover the 
lowest order perturbative anomalous dimension of the field in four dimensions has also been obtained jjjj from the 
WH equation by making use of a sharp cut-off to integrate out the ultraviolet modes. 

In this paper we critically reconsider the set of coupled equations for Uk(<i)) and Zk(4>) and more generally discuss 
the derivative expansion in the framework of the WH equation. In particular, after a brief review in Sect. II of the 
WH method and the derivation of the LPA, in Sect. Ill motivated by the result of we check the reliability of the 
equations derived in that paper and evaluate the anomalous dimension of the <fi 3 theory in D — 6 dimensions to the 
order 0(<?§) and, afterwards, the anomalous dimension of the scalar theory at the Wilson-Fisher (WF) fixed point in 
D < 4 dimensions. 

In Sect. IV we move to a critical revision of the derivation of this system of coupled differential equations for Uf. and 
Zk ■ This gives us the possibility of reconsidering the longstanding problem which affects the implementation of the 
derivative expansion of the WH equation. As soon as a non-constant field is considered, the differential equation for 
Zk(4>) (and, more generally, the equations for the higher order coefficients in the derivative expansion), is affected by 
the presence of some non-analytical terms that apparently spoil the very differential nature of the equation. The origin 
of these terms is analysed in detail and the derivation of the equation for Zk in the momentum space is presented. 
We shall see that under a specific assumption of smoothness of the background field, it is possible to derive the 
differential equation and that the dangerous terms are actually negligible. The fate of these non-analytical terms 
when the sharp cut-off is replaced with a smooth one is also discussed and we shall show that for a sufficiently slow 
fluctuating background field, both sharp and smooth cut-off produce the same equations. In addition we show that 
the assuption about the relative weight between kinetic and potential terms in the action is crucial in the derivation 
itself of the Zk equation. Our conclusions are summarized in Sect.V. 



Let us call Sa[3>] the action of a scalar field $ which contains all the Fourier components <fi q such that < \q\ < A. 
For the sake of simplicity in this paper we shall only consider the single component scalar theory.Usually we want to 
describe a physical process as the scattering of two bosons with momenta p^ such that p << A. As already mentioned 
the origin of the difficulties in perturbation theory is due the fact that we have to take into account the contributions 
to this process coming from all the Fourier modes between p and A. Actually the action 5a [^J, often called the 
wilsonian effective action at the scale A, is more appropriate to describe the physics of the system at scales p ~ A. 
At a scale p << A it would seem more convenient to directly deal with an effective action S p [<j>] that describes the 
physical phenomena taking into account only those scales around p. The Renormalization Group Method provides 
this effective action. Let us define S p [<j)] , in the Euclidean version of the theory, through the equation 



where we have split the original field $(x) into a background field 4>{x) containing only the modes between zero and 
p and a fluctuation field (,{x) containing those between p and A, 



S p [</)], called the wilsonian effective action at the scale p, contains, within its parameters, the effect of the interactions 
among the modes in the range [p, A] and the modes below p and therefore it is the effective action we were looking 
for. The difficult task here is to perform the integration over the high frequency modes ^. The above integral can 



II. LPA AND LOOP EXPANSION 




(1) 



$(x) = <j)(x) +CO). 



(2) 
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be approximated by means of the loop expansion which amounts to a series expansion in powers of fi. The one loop 
approximation gives the O(h) correction to the tree level result S p [<f>] = S&[cf>] and corresponds to a gaussian integral 
over the fluctuation £ obtained expanding S\ [<f> + C] U P to C 2 0] ■ 

We may now pose the question in a different way and ask how does the wilsonian effective action Sk evolve once 
we integrate the modes in the shell [k — Sk, k] with an infinitesimal 5k to get S k ~8k- In this case the gaussian 
approximation, i.e. the one loop result, becomes exact J2j. The problem of the evaluation of the effective action at 
the scale k in Eq. (|l|) can then be turned into the problem of solving a first order differential equation for Sk [4 1 } with 
respect to the variable k. The action Sa[$] is the the boundary condition for this equation. This is the Wegner- 
Houghton method. Actually the problem of solving this differential equation for Sk is as difficult as to perform the 
integral in Eq. (Q) and we have to resort sooner or later to some approximation if we want to make any progress. As 
a first step let us go back to Eq. (1) which, replacing A with k and p with k — Sk, is the defining equation for SkSk 
in terms of Sk 



[D(]e 



(3) 



Our second step is to insert in the above equation an ansatz for Sk • At this point we make the assumption that 
a good approximation to Sfc[$] is provided by the derivative expansion (D is the number of dimensions) 



d D x 



(4) 



By inserting Eq. (Q) into (||) we should be able to derive an infinite system of coupled differential equations for the 
coefficients functions Uk , Zk , Yk , ■ • • • 

The lowest order approximation in the derivative expansion is the LPA. Replacing in Eq. (^) Zk = 1, Yk = 0, • • •, 
and then considering a constant background field <f>(x) = 4>q we get from Eq. (||) an evolution equation for Uk{4>o) 
only 



d 

k—Uk(M 



hk D N D , 



ln(k z + U k (<p ))- 



(5) 



Here the ' means derivative with respect to 0o and the result of the angular integration is Njj — 2/ 1 {{Ait) d / 2 T{D /2)). 

This equation has been found again and again (see for instance Jl4|,p[) and the consensus on this equation is 
unanimous. 

The equation above is a non perturbative evolution equation for Uk- Suppose that the potential Uk{<po) has a 
polynomial expansion and for the sake of simplicity let us require the Z(2) symmetry 0o — > — 0o 



(6) 



By taking 2n times (n — 0, 1, 2, • • • ) the derivatives of Eq. with respect to <po at 4>q = we get an infinite system 
of coupled equations for the coupling constants g2n(k): 



k§- k g 6 (k) 



Tik D N D g A (k) 



hk D N D 



9l(k) 



2 +g2(k) 

ge(fc) 
k 2 +g 2 (k) 

fls(fc) 

k 2 +g2(k) ^(k 2 + g 2 (k)) 
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(k 2 + g 2 (k)) 2 
g4(k)g 6 (k) 



30 



{k 2 + g 2 {k)f 



(7) 



A diagrammatic interpretation of the above equations is straightforward, provided one identifies g 2n with the 2n 
external legs vertices and the denominators (fc 2 +g 2 {k)) m with m propagators joining the various vertices; in fact each 
r.h.s.in Eqs. (Q) represents the sum of all the one loop diagrams with fixed number of external legs that can be arranged 
combining the g 2n . It is then clear that the above system is an approximation to the infinite set of Schwinger-Dyson 
equations for the Green's functions at zero external momenta. Of course as it is true for the complete Schwinger-Dyson 
system, Eq. (Q) is indeterminate. To make contact with perturbation theory, more precisely with the h expansion, we 
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can seek for solutions where each coupling constant is developed in an S-power series gi n = g^J + hg^J + h 2 g^J + • ■ ■. 
As it is well known in this case the solution of the Schwinger-Dyson equations and a fortiori of our system becomes 
unique and it may be easily verified that the ultraviolet behavior of the coupling constants flow at 0{Ti) coincides 
with the usual one- loop result. Let us consider for instance the one loop contribution to the perturbative (3- function 
for the </> 4 theory. It is recovered in the ultraviolet regime, k 2 » gi, by retaining in the r.h.s. the 0(h ) values of 
the couplings, that is their boundary values at k = A, .94(A) = <7 4 , (76(A) = gs(A) — ■ ■ ■ = 0. However Eqs. (]?]) show 
that <76 is 0(174) and then g$ in the r.h.s. of the /3-function of 54 provides a 0(lj\) effect. It is then clear that higher 
loops perturbative contributions to the /3-function are present in Eqs. (Q) due to the fact that all vertices in the r.h.s. 
are fc-dependent. It is also easy to realize that, in order to recover the full two-loop 174 /3-function, it is necessary to 
go beyond the LPA since the ge and gi contributions are not sufficient to get the complete 0(gf) effect. 

Going back to Eqs. (]?]), in order to better exploit the non-perturbative character of this system, a different kind 
of truncation is needed. A first step toward its non-perturbative analysis has been taken in Jig ]. However, due to 
the appearance of non-trivial saddle points in the renormalization group equation studied in that paper, the results 
should be taken quite cautiously. 

The first step beyond the LPA is discussed in the next section. 



III. DERIVATIVE EXPANSION AND ANOMALOUS DIMENSION 



By turning on the scale and field dependence in Z k we allow for a non-trivial lowest order derivative term in the 
action (Q), thus obtaining the first improvement to the LPA where it was set Z = 1. In this case we have to deal 
with two coupled equations for Uk and Zk which can again be reduced to an infinite set of coupled equations if one 
assumes a polynomial expansion for both U k and Zk- 

In [ fL^ ] a specific procedure has been carried out in order to determine a differential equation for Zk in the sharp 
cut-off limit and it has been shown that in D = 4 for the O(N) symmetric tfi 4 theory, the perturbative anomalous 
dimension of the field to 0(g~ 2 ) can be obtained from that equation. However it was also noticed that, following 
a different and, in principle, correct procedure, a different equation for Zk is obtained. In practice the gaussian 
integration of the fluctuation ((x) in Eq. (|J) yields a contribution to Sk-Sk proportional to ln[(S 2 Sk/S^5^)\^ =( j > ] 
and both procedures require to expand the logarithm in order to derive the differential equation for Zk- In one case, 
following the steps outlined in ||l6| for the one component theory, we split the logarithm (again ' indicates derivation 
with respect to the field) 

ln[(5 2 S k /5$5$)\^} = ln({Z k ) d^ + (C4)o + a) (8) 

where (Z k )o and {Uk) are evaluated at the constant field configuration (fto and A is the x dependent part of the 
propagator; the expansion is performed taking A small with respect to the non-fluctuating part. In the other case, 
according to the procedure developed in |l7j , 

ln{(6 2 Sk/6<S>6<S>)\^} = ln(z k k 2 + u' k ' + Z k (d^ - k 2 ) + c) (9) 

where Zkk 2 is added and subtracted and the remaining part of the propagator C is proportional to derivatives of the 
field. In the latter case the expansion is made requiring that C and the difference Zkid^d^ — k 2 ) are small. This is 
justified due to the field derivative term in C, which is small within the derivative expansion framework, and to the 
fact that d fl d fJ ,6(x — y) yields in the infinitesimal integration shell [k — 5k, k) a factor k 2 which makes the difference in 
Eq. (^) vanishing. It is not surprising that two expansions obtained for different choices of the expansion parameter, 
can be not straightforwardly comparable. As a consequence of the different nature of the expansions in Eqs. (|J) and 
(^J) , we end up with two different equations for Zk . We will come back to this point in the next Section. 

The point of view chosen in was to test the reliability of these expansions by comparing the values of the 
anomalous dimension of the field 77 at the lowest non-vanishing perturbative order obtained from the equations with 
the result of the usual perturbative method. It turned out that the equation obtained from the expansion @ provides 
a value of 77 in agreement with perturbation theory in D = 4. 

In order to support the result of Jl^] and to show that it is not an accidental coincidence, one can easily repeat the 
calculation of the anomalous dimension performed in |i~2| for the case of the cubic theory 4> 3 in D = 6. Obviously eqs. 
([|[7]) are not consistent with the symmetry of the latter theory, but the complete equation for the potential (||) and the 
equation for Zk, deduced in p2| , are still valid because they do not require any assumption on the internal symmetry 
of Uk and Z k - The explicit form of the Zk equation as deduced from the expansion of Eq. (j^) in D dimensions is 
(A = Z k k 2 + U' k ' and % = 1) 
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d_ k D N D( zl 2Z' k A' Z' k \ 2 2Z k A' 2 8Z' k Z k A'k 2 2Z 2 A' 2 k 2 ^ 

dk k ~ 2 \A A 2 DA* + 3A3 + 3DA3 DA 4 J [ ' 

and, in order to get a perturbative estimate of the anomalous dimension, we replace, in the r.h.s. of Eq. (flO|), U k and 
Z k with the corresponding bare quantities, namely 7? 3 </> 3 /3! and 1 (we perform our calculation in the ultraviolet regime 
where the mass term in the potential can be neglected). The integration of Eq. ([l(]), from A down to a generic value 
k < A, is then straightforward yielding, in D — 6, Z k — ~g 2 / (3847r 3 )/n(A/fc). The anomalous dimension is immediately 
deduced from this ^-independent expression of Z k 

(11) 

and it is in agreement with the perturbative diagrammatic computation |Tg| ] . Actually the usual perturbative approach 
in this case is very simple since it involves a one-loop computation and, correspondingly, we need only one integration 
of Eq. ([To]) to get the answer, differently from the D = 4 case where, on one side a two-loop diagram is to be evaluated 
and, on the other side, a two step integration of Eq. ([!(]) is required to determine the lowest (non-vanishing) 
perturbative contribution to 77. However, the above calculation is similar to the one performed for the <fi theory in 
four dimensions since in both cases the anomalous dimension is perturbatively expressed as expansion in powers of 
the marginal dimensionless coupling appearing in the potential. The running of these couplings with the scale k is 
only logarithmic thus justifying their replacement with the bare constant g 3 for D — 6 and g 4 for D = 4. The same 
argument holds for the field independent part of Z k which is dimensionless and can be replaced with its bare value 1. 

The further important step is to check whether this equation is able to reproduce the value of rj at the non-gaussian 
(WF) fixed point which appears in the scalar theory below four dimensions. 

In order to get informations at a fixed point it is convenient to express all dimensionful quantities in terms of the 
running scale k and rewrite the coupled equations (|5|,[To|) in terms of dimensionless variables which can be conveniently 
introduced through the relations t = ln(k/A), x = k^- D -^/ 2 {2/N D ) 1 / 2 (j), u(x,t) = 2k-°U k /N D , z(x,t) = WZ k . In 
these relations powers of Nd/2 have also been included in order to get a simpler form of the differential equations; 
indeed the constant a — Nd/2 disappears from Eqs. (|E|),(|l0|) after the replacement U k — > all k , <fi — > \/a<p, Z k — > Z k . 
Furthermore, since we are interested in carrying out a numerical analysis of the problem, instead of directly attacking 
Eq. (||) which contains a logarithm, we shall consider the corresponding equation for the derivative of the potential, 
as already performed in 0. Actually the presence of the logarithm causes a stiffness of our set of equations which 
is a source of many numerical drawbacks. Therefore, by rearranging the equations in terms of the derivative of the 
scaled potential f(x, t) — du{x, t)/dx we finally get the two coupled differential equations (a = z + f and this time ' 
means derivation with respect to x) 

v__ D/ . + £z|±2>c, +l/0 _*±£ (12) 

at 2 a 



dz (D-2 + rj) /z" 2z'a! z' 2 2za' 2 8z'za' 2z 2 a' 2 \ 



and the fixed points correspond to t- independent (df jdt — dz/dt — 0) solutions of Eqs. ( |12| , |13D . In the following 
part of this Section we shall consider only t- independent solutions of Eqs. ( )l^ , |l3| ). 

As discussed in [p|j8|JT9[| , Eq. ( |l2| ) alone allows to determine the fixed point structure of the theory and specifically the 
appearance of the WF fixed point is shown for D < 4. It is easy to check that the gaussian fixed point, corresponding 
to / — 0, z — const ^ 0, 7] — 0, is a solution of ( |l2| , [i"3| ) for generic D, whereas the determination of the non-gaussian 
fixed point with the corresponding value of rj requires a numerical analysis. 

Eqs. (1203) are solved by requiring the usual normalization of the kinetic term in the action and two constraints 



on /' and Zctt x — which preserve the Z(2) symmetry 

*(0) = 1 /'(0)=0 z'(0) = (14) 



As explained in jS|, from eqs. ( |12| , |13| ), one can easily determine the asymptotic behavior of the solutions f(x) and 
z(x) for large x, up to two (one for each solution) multiplicative constants, c/,c z . Thus we have enough boundary 
conditions to integrate our equations and to deterimine the three unknown constants Cf,c z and rj. In order to deal 
with a two point boundary problem, we have used, as also suggested in [^j, the 'shooting method' embedded in a 
Newton- Raphson algorithm po[ . The results shown below have been obtained by requiring that the difference between 
the constraints in (|14| ) and the shooting variables is less than 6 = 10~ 7 . 
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In D = 3 the value obtained for the anomalous dimension is 77 — —0.071, quite different from the world best 
determination 77 = 0.035, quoted in B, or the value determined through the e-expansion to the order 0(e 3 ), 77 = 0.037 
(see pl|). The result is very poor even if compared to the one obtained starting from the equation for Zk derived 
resorting to a smooth cut-off fa]. 

However it should be noticed that the situation is less bad if D increases. This is illustrated by comparing the value 
obtained for 77 at the WF fixed point in D = 3.4, 77 = 9.74 10 -3 with the estimate of 77 obtained from the e-expansion 
g§| to 0(e 3 ), r] = 10.70 10~ 3 , and <3(e 4 ), 77 = 9.62 10~ 3 . In D = 3.6 we get 77 = 4.30 10" 3 to be compared with 
the value from the e-expansion to 0(e 3 ), 77 = 4.16 10 -3 and to <3(e 4 ), 7; = 3.95 10~ 3 . In D = 3.4 and D = 3.6 the 
differences between the various estimates are below 10%. In Figs. 1,2,3 we show f{x) and z(x) obtained at the WF 
fixed point for D — 3, 3.4, 3.6. 

All curves for f(x) in Fig. 1 intersect the x-axis at non-vanishing values of x, which correspond to non-zero minima 
in the potentials; the curves z{x) in D = 3.4 and D = 3.6 are plotted in Fig. 2; analogously to what is found in || 
for the corresponding variable in the smooth cut-off framework in D — 3, z(x) has a maximum and then decreases 
for large x (even the dashed curve, although less evidently, smoothly decreases after a maximum at about x = 8). 
Conversely, as shown in Fig. 3, z(x) in D = 3 is increasing: no maximum has been found even enlarging the x range 
to the limits allowed by the integration routine (in practice the upper limit can be pushed up to about x — 12). 
Moreover, close to x — 0, z(x) < 1 and the curve has a minimum which is not present in the curves in Fig. 2. 

In order to justify this peculiar behavior, we notice that the leading terms within the brackets in the r.h.s. of Eq. 
( |l3| ) are the ones proportional to z and not z', namely (2za' 2 )/(3a 3 ) — (2z 2 a' 2 )/(Da 4 ), that is, as long as D is not 
close to 3, these two terms are dominant since z » z' but in D — 3 and in the ultraviolet limit, when a ~ z ~ 1, 
they practically cancel out and the behavior of Eq. ([[3]) is sensibly modified. 

In conclusion, only close to three dimensions, where the infrared effects become more and more important, Eqs. 
(|l2],[l3]) fail to reproduce the known results about the anomalous dimension. 



IV. WH EQUATION AND NON-ANALYTICAL TERMS 



We now move to another point which, in some sense, is preliminary to the analysis of the coupled equations for the 
various parameters entering the local action, namely the rise of some undesirable terms as soon as one goes beyond 
the LPA in the derivative expansion. The presence of non- analytical terms has been recognized since long time pH in 
relation with the use of a sharp cut-off and recently reconsidered in |0,^,^3). Here, in order to get a clearer insight 
into this problem, we shall study the origin of these terms working out the various steps for the determination of the 
differential equation for Zk- 

Let us go back to Eq. (||) and (Q). The LPA approximation corresponds to setting Zk — 1, dropping Y% and all the 
others higher order derivative coefficients and finally restricting to a constant background <j) = </>o . In order to get the 
equation for Zf., we have to release the condition Zk = 1 and collect all terms proportional to d^ifid^ in the r.h.s. of 
Eq. (||) and, to this purpose, it is necessary to retain a non-constant background cf>(x) — <fio + (p(x). Let us choose 
for the non-uniform component of the background (p{x) a single mode with |g| << k (the reason for this choice will 
become apparent later) 

i P {x) = ^={ ( p g e i ^ + ^ q e-^ x } (15) 

Here V is the volume factor. As before, the fluctuation C(^)i which must be integrated out, contains Fourier modes 
only in the shell [fc — 8k, k] 

CO*) = -jrrzZ^^ ( 16 ) 

The square bracket in Eq. ( |l6| ) , [p] , is a reminder that the sum is restricted only to those values of p such that 
k - 5k < \p\ < k. 

Therefore the original field is split into three parts, $(x) = 4>q + f(x) + C( x )- ln order to pick up contributions up to 
the second derivative term ^$9^$, the functions Uk{4>o + (p(x) + ((x)) and Zk{4>o + f{x) + ((x)) must be expanded 
around </>o up to ip 2 {x), and since the gaussian functional integration is exact, as mentioned before, we only need to 
consider terms in ((x) up to 0(£ 2 (x)). Due to the condition q << k the linear terms in £(x) drop out after the spatial 
integration. Let us consider now the terms proportional to C 2 { x )- To illustrate the various steps made to derive the 
equation and the related problems we do not need to consider here all of them. It will be sufficient to consider for 
instance those coming from the expansion of the potential. Performing the spatial integration with the help of Eqs. 
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|15|,P-6|) (in the following we shall consider, with no loss of generality the four dimensional case D — 4 and, in order to 



avoid any misunderstanding of notations, the n times derived potential is indicated here as C/i"') we have 



d 4 x( 2 (x) 



U, 



(2), 



(17) 



u ^ W)c 2 (*) = ^ o) 



2VV 



Vq E CpC-p-g + f-q E CpC-p+q 

v [p] [-p-q] [p] [-p+q] ' 



(18) 



(4), 



AV 



d 4 xip 2 (x)( 2 (x) 

Pq E CpC-p-2g 
[p] [-P~2g] 



2 VqV-q2^ 

[P] 



CpC— p 



E CpC -p+2g 

[p] [-p+2g] 



(19) 



One more comment about the notation. All sums above are single summations over p and the range spanned by p 
in each sum is obtained by requiring that all quantities in square brackets are constrained into the shell [k — Sk,k]. 
For instance in the first sum in Eq. (fl8|), [p] [— p — q] indicates the double constraint (k — 5k) < \p\, \ — p — q\ < k, 
which in turn determines the values of p selected in the summation. Then, as soon as we turn on a non-uniform 
background f(x), we end up with a deformation of the integration region. The new integration region for p is now 
given by the overlap of two shells with one of the two depending on q. Only at q — the two shells are both reduced 
to the original one. The reason for the appearance of the non-analitycal terms has to be traced back to this boundary 
effect. 

In order to deal with shorter expressions we limit ourselves to the case of a (^-independent Z k . Then, with the help 
of Eqs. ( ^7| , ^8|Jl9| ) , we can easily expand the action Sk[4>o + f + C] around </>o and perform the quadratic integration 
in Q obtaining (fi = 1) 

exp{-S k -5k[4>o + f]} = exp{-Sk[4> Q + </?]} exp < ZnG _1 (p) > 

W J 



x{l + %W, 



Y, G(p)G(p + q)+ J2 G(p)G(p-q) 
[p] [p+q] [p] [p-q] 



(20) 



where we have introduced the propagator-like notation G 1 {j>) = Zkp 2 + U k ((f>o) + l/j^' ((/)o)(p q tp- q /V . 

Since we are just interested in the Z k evolution, we again neglect cubic and higher powers of tp and expand G(p), 
G(p + q), G(p - q) around q = up to 0{q 2 ), obtaining (A(p) = Z k k 2 + U { k 2) (cf) ) ) 



Sk-5k[<l>o + ¥>] = Sfcl^o + <p] + ^E ln ( A (p)) + 



E 

[p] [p-q] 



-PqP-q* 



2v w-^m 

[p] 



E 



i 



i - 



q 2 Z k A(p-q) 2 Z 2 



A 2 (p) L A(p) 



V[p] [p+q] 

q 2 Z k A(p-q) 2 Z 2 



A 2 (p) 



A 2 ( P ) L A{p) 



A 2 ( P ) 



(21) 



At this point if we just neglect the boundary problem and everywhere integrate the momentum p within the shell 
[k — 5k, k], from the terms in Eq. ( pl| ) not depending on tp we get 



d k 4 
k—Uk^o) = -—ln{Z k k 2 + U' k \4> )) 



dk 



(22) 
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This equation replaces Eg . (^|) in D — 4. From the terms proportional to (p q cp- q , we get the second derivative with 
respect to 0o of Eq. (§2f. This is simply because those terms come from an expansion of Eq. (^2|) evaluated at 
4>{x) — 4>q + ip(x) around 4>a- 

Finally if we collect the terms proportional to q 2 ip g Lp- q we get the evolution equation for Z k 



a fc 4 l z k u£\^) zluj 3 \M 2 P 

dk k ~ 16tt 2 I A{kf A(k) 4 ) 1 ' 



In the absence of the boundary problem and under the assumption of a field independent Zk , (|22| ) and fl23|) are the 
system of coupled differential equations for U k and Zk- In the mo re g eneral case of a field dependent Zk by following 



the same steps as before we find that the equation that replaces (23) is 



k jL 7 fc 4 /< 2z' k A> z' k fc 2 z k A' 2 z' k z k A>e Z 2 A' 2 k 2 ^ 

dk k 16tt 2 \ A A 2 AA 2 A* A 3 A 4 J [ ' 

Several comments are in order. This equation has already been derived by following different steps in (UJ2^]. As 
already noticed in |t2| it is substantially different from Eq.(jl^), namely both equations contain the same kind of terms 
but with different numerical coefficients. In the previous Section we have already explained that to derive Eq. (|k|), 
an expansion where the "kinetic" term Zk{df_ L d^ ~ k 2 ) is considered small w.r. to the "potential" term Zkk 2 + U k is 
used. The expansion that leads to Eq.(^4j) is of a completely different nature. In this latter case the potential term is 
considered small w.r. to the kinetic term. It is then not surprising that we get different results. It is well known that 
no one of these expansions can be considered as definitely superior to the other. The choice of one of them is rather 
dictated by the physical problem under investigation. 

It may be stressed at this point that the fact of having derived two evolution equations for Z k , namely Eqs. (^4|) and 
( |l0| ) , is not the consequence of having used two different definitions of this parameter (which is actually introduced 
by the derivative expansion (|J) of the action whose evolution is determined by the integration in Eq. (||) , with the 
boundary condition Z\ = 1). Rather, it is the consequence of having used two different approximations, which in 
turn are necessary since the full problem cannot be solved exactly. More specifically, as we have just seen, we need 
to insert explicitly a non-uniform ( and slowly varying, so that higher derivatives terms in the derivative expansion 
can be safely neglected) background field to read the differential evolution equation for the coefficient Zk- Now the 
fluctuation operator Sip(xf84>(y) ^ s obviously no longer diagonal in momentum space. To determine then a differential 
equation for Zk , we need a further expansion in which a piece of the propagator is considered small with respect to the 
rest of the propagator. In the derivation of Eqs.(|i"o[) and ( |24| ) above this additional expansions correspond respectively 
to the physical case where the kinetic energy is small with respect to the potential energy and viceversa. 

The method exposed above to obtain Eq. ([!(]) is different from the methods used respectively in and [p5| . 
Nevertheless the expansion adopted in the three cases, namely a small coupling expansion, is the same and this 
explains why these methods all lead to the same equation. 

Finally Eq. (^4|) is obtained by neglecting the boundary effects due to the deformation of the integration shell. Let 
us now go back to Eq. (|2l] ) and work out explicitly the consequences of this boundary distortion. To illustrate this 
problem we just need to consider one term. Let us take for example the term 

Uk (3) (M 2 \- l Z k q 2 

[p] [p+q] y ' v ' 

To further simplify the point, let us write down explicitly the above indicated sum in D = 1 dimensions. The only 
additional complications that would appear when considering higher dimensions are some factors coming from the 
angular integrations. By trivially replacing the sums with the integrals we get as a result of the overlap of the two 
shells (L is the volume of the system) 



1 Z k q 2 _ a L / f- k+5k dp , f K dp 

-= Z k q — < 



y 

MtV 2(p) A{p) (z kP 2 + u { k 2) m 3 J^+i(z k p 2 + u { k 2 \<P)) 

= Z k q 2 ±- 2 — Jsk-q) (26) 



(z k k 2 + u^>m 



It must be noted that the companion term in eq.(|21j), which differs from the one in Eq. ( |2q ) only for the integration 
limits, gives the same result and we do not get any cancellation of the term proportional to q. The above result 
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explains our previous statement that the appearance of non analytic terms in q 2 is simply due to the distortion of the 
integration domain induced by the non uniform background configuration. 

In Eq. ( p6| ) a term proportional to q appears which is not proportional to 5k. This in turn seems to imply that 
the very possibility to establish a differential equation for Z k is in danger and actually the presence of such terms 
has been considered as a serious drawback that blocks the way to any application of the WH method beyond the 
LPA approximation. That consideration has convinced some authors that to overcome this problem a smooth cut-off 
rather than a sharp cut-off should be used. We shall come back to the smooth cut-off later. 

The above derivation actually teaches us that there is no serious drawback as far as we correctly interpret our 
equations. The undesired non-analytic terms that appear on the r.h.s. of the defining equation for Sk-Sk and that 
are not contained in the original derivative expansion ansatz, always appear in combination with 5k as in the example 
above. This means that as long as we keep q << 5k, that is as long as the non-uniformity of the background field 
is small compared to the width of the shell within which the modes arc treated as independent, these terms can 
be coherently neglected and the equations above are the correct WH RG equations in this approximation. This is 
actually one important result of our analysis. The very possibility to implement the RG trasformations as a system 
of differential equations for Uk and Z k ( i.e. to consider such a truncation to the infinite system of equations for all 
the coefficient functions that appear in the derivative expansion as a good approximation) is intimately related to the 
nature of the background field <p(x) that we consider. It has to be sufficiently smooth over the scale 5k otherwise the 
differential equations simply cease to be valid. The important point here is that not only we learn that it is possible 
and perfectly legitimate to establish differential equations that implement the infinitesimal RG transformation even 
beyond the LPA approximation, but we also obtain the limit of validity of these equations. 

It is worth at this point to pause for a moment and review all the steps and approximations involved in the above 
derivation. First note that the wilsonian renormalization group transformation gives: 

S k -sk(<l>) = S k sk(cl>) + 5k(- ■ ■) + 0((5k) 2 ) (27) 

From the derivation of Eq.(|2^) it is clear that 5k is the momentum range above which the modes are considered as 
independent, that is any momentum function f(p) that appears in the above equation is practically constant within 
this range. This in turn means that we can neglect the 0{{5k) 2 ) terms or, in other words, that the above equation 
can be rewritten as a differential equation. This is what is meant when the above equation is referred as an "exact 
differential renormalization group equation" and this is also the meaning of the equivalent statement saying that the 
gaussian integration is "exact". Actually all that means that the 0{{5k) 2 ) terms can be consistently neglected. 
Next we approximate S k by a derivative expansion. We have seen that as soon as we go beyond the lowest order and 
allow for a non constant background field we encounter singularities that make the formal mathematical limit 5k — > 
ill defined. We have seen above that this is none of a problem as far as the scale at which singular terms appear is 
kept far from the "resolution scale" 5k. This is nothing but the condition q << 5k. 

The two points above can be easily understood by considering the following hydrodynamical example. To give 
a differential form to the equation for, say, the density p of a fluid we need two conditions to be satisfied. The 
infinitesimal volume d 3 x has to be "sufficiently small" so that the density as well as any other macroscopic quantity 
can be considered as constant within this volume. At the same time it has to be "sufficiently large" so that a 
macroscopic number of molecules is contained within this volume. Under the above conditions we can write the 
evolution equation for the density in a differential form because the 0((d 3 x) 2 ) contributions can be neglected and 
because due to the second condition we can neglect the singularities coming from the molecular scale a , once 
our "resolution scale" dx has been taken much bigger than a : dx >> a. The complete analogy between the 
approximations involved in the derivation of our equations and the hydrodynamical example should be, by now, clear. 
Obviously the differential equation for p obtained this way becomes less reliable and ultimately totally wrong when 
we want to describe physical phenomena whose resolution scale approaches the molecular scale a. The same warning 
applies to our case when the background field is not sufficiently smooth within 5k, and this happens when q — > 5k. 

Let us now turn our attention to the smooth cut-off procedure. To illustrate the point we can again reconsider 
the previous example and implement the constraint in Eq. (|2^) by means of differences between theta functions in the 
following way 

\ - 1 Z k q 2 2 L [+°° e (p 2 , k 2 , (k ~ 5k) 2 )Q Q ((p + q) 2 , k 2 , (k - 5k) 2 ) 

MW y ' KJ {Z kP 2 + U k ' (0o )) 

where we have defined Qo(p 2 ,k 2 , (k — 5k) 2 ) = 6(p 2 — (k — 5k) 2 ) — 9(p 2 — k 2 ). 

The effect of a smooth cut-off is implemented here by replacing the Oo(p 2 , k 2 , (k — Sk) 2 ) function with a smoothened 
version Q e (p 2 , k 2 , (k — 5k) 2 ) where e has the dimension of a momentum and is a new cut-off that we introduce in 
the theory; in the limit e — > it should give Qq- It is easy to see that all the non-analytic terms disappear. In fact 
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by expanding the function c ((p + q) 2 , k 2 , (k — 6k) 2 ) around p 2 we get (from now on we omit k and k — 5k in the 
argument of the O e function) 

e e ((p + q ?) = e.b 2 ) + (2 Pq + q 2 ) + I^jp (2w + <zf + • • • (29) 

Being now the integration region symmetric all contributions from odd powers of p vanish and consequentely odd 
powers of q are absent. As we shall see in a moment we have gained nothing from having made the non-analytic 
terms disappear. Let us focus on the function e and regard it as a function of \p\. The scale e gives the size of the 
region over which O e as a function of \p\ changes significantly from zero to one. This means that its first derivative 
with respect to \p\ is 0{\) in the two regions around k and k — 5k, and zero everywhere else. Analogously its second 
derivative is O(t-) in the same region and zero everywhere else, and so on. 

We could expect that e is introduced just as an intermediate step and that the final results are obtained after sending 
e to zero and they arc finite. The above result shows that this is not the case. As a consequence of the introduction of 
the smoothening cut-off e we have generated i divergences which have already been noticed in |t],^2|,^3| . This could 
appear at a first sight as a disaster and again we might wonder whether it is possible to establish differential equations 
for Uk and Zk ( and more generally for the coefficient functions of the derivative expansion) under these conditions. 
Actually we can repeat here the same kind of considerations that we have done in connection with the appearance of 
the non analytic terms. After performing the integration in the momentum p we see that 

i) the terms proportional £ (p 2 ) give contributions 0(5k). This is because the integration region where the integrand 
is significantly different from zero is the shell [k — 5k, k]; 

ii) the terms proportional to - give contributions proportional to (ijli (f) Qi an d so on. This is because 
the region over which the integrand is significantly different from zero has now a width of e and we get then a factor 
e coming from the width and a factor - coming from the derivative; 

iii) the terms proportional to - give contributions proportional to ff^QS f§) (<0^' (^) (<0^ ! (^) (*)^' 
and so on. This is because the region over which the integrand is significantly different from zero has a width e as 
before but this time we get a contribution O(t-) from the second derivative of C . 

Let us now comment on the above results. First we make the choice e ~ q, that is a very good choice for e 
once we remember its role of smoothening parameter around the points k and k — 5k and that we want to integrate 
essentially only the modes within the shell. It is now apparent that as long as we keep the scale q « 5k the dominant 
contributions come from the terms proportional to e (p 2 ) and all the contributions coming from the derivatives of 
9 c (p 2 ) can be coherently neglected and the final output is practically equivalent to the sharp cut-off one. 

An important comment has to be made at this point. The particular choice of the smoothening function e that we 
have done above, more precisely the fact that we have chosen it to be an even function of \p\, is responsible for getting 
only even powers of q in the above results, in other words we do not get non analytic terms in q 2 . But this is absolutely 
irrelevant because what matters to establish the differential equations is the condition q « 5k irrespectively of the 
fact that we have odd or even powers of q. In fact had we chosen G c to have an odd dependence on \p\ we would have 
also found odd powers of q. Nevertheless under the conditions q << 6k and q << k, we can equally well neglect all 
the terms apart from the first and again recover the differential equations. 

What we have just seen proves that actually there is no conflict between the sharp and the smooth cut-off approach. 
These two procedures, once correctly implemented, give precisely the same results as it should be expected. We have 
also learned under which conditions the couple of differential equations for Uk and Zk are valid. 

Another additional comment peculiar to the smooth cut-off procedure has to be made. Due to the disappearance 
of the non-analytic terms there has been in recent years some preference for the smooth cut-off implementation of 
the differential RG transformations w.r. to the sharp cut-off. As we can easily see from the points (i) (ii) and (iii) 
a large value of e looks even more efficient in suppressing the undesired terms and we could be led to the conclusion 
that we can accomodate such large values of e within the smooth cut-off approach. Some results have been recently 
derived by making use of this apparent better flexibility of the smooth cut-off versus the sharp cut-off. Actually the 
condition of having a small e ( where small means not too big compared to 5k) is necessary to be coherent with the 
very strategy of the RG method as exposed before. Substantially only the modes within the shell have to be integrated 
out. By considering larger values of e we move toward the independent mode approximation. The conclusion is that 
those results that have been obtained within the framework of the smooth cut-off procedure with the help of cut-off 
functions whose typical width is of 0(k), where k is the UV cut-off, have to be taken with a grain of salt and the 
improvement w.r. to the perturbative result (independent mode approximation) is not very much under control. 
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V. SUMMARY AND CONCLUSIONS 



We have carefully analysed the derivative expansion for the effective action and shown how the exact rcnormalization 
group equations for the coefficient functions Uk and Zk are obtained. The most important lesson we have learned is 
that the system of coupled differential equations for them can actually be established provided the background field 
around which the quantum fluctuations are integrated is sufficiently smooth. The width Sk in Fourier space measures 
the momentum range within which the modes are treated as independent, meaning that any function f(p) that has 
to be integrated in the shell [k — Sk, k] is considered to be constant within it. The key point of the method is that this 
shell has to be on the one hand sufficiently small, i.e. Sk << k, where k is the UV cut-off, so that the feed-back of the 
higher energy modes on the lower ones is correctly taken into account. On the other hand, and this is the crucial point 
that we want to emphasize here, it has to be sufficiently large in such a way that within this range the background 
field can be considered practically flat, i.e. q « Sk. Having made clear this point we have seen that the non-analytic 
terms are no more source of problems and the differential equations can be safely established under these conditions. 
We have also proven that the introduction of a smooth cut-off practically does not change our conclusions. In fact, 
under the conditions quoted above, sharp and smooth cut-off produce precisely the same results. A word of caution 
has to be said regarding the results that have been obtained within the smooth cut-off procedure by allowing the 
smoothening scale e to take values O(k). In this case the danger is that the width of the shell Sk gets infinitesimal 
w.r. to the smoothening scale which means that the modes are practically all treated as independent and then there 
is no control on the possible improvement on the perturbative results. 

Another important point we have learned from the above analysis is that depending on the physical problem at 
hand different expansions can be envisaged and this lead to different equations. The typical situations illustrated are 
the one in which the kinetic energy term is considered small compared to the potential and the opposite one. These 
two cases respectively produce Eq.(|lO|) and Eq.(|24|) above. 

One third important result of our analysis is that ( see Sec. Ill) the study of the coupled system of equations for 
Uk and Zk has shown that we are unable to reproduce the value of r\ at the Wilson-Fisher fixed point in D — 3 
dimensions. Our previous results teach us that the reason of this failure is not to be searched in an intrinsic weakness 
of the sharp cut-off procedure versus the smooth cut-off. Presumably the reason is that we enter a region where the 
equations themselves are no longer valid. An indication for that can be seen in the behaviour of Zk- When it deviates 
from Zk = 1, then less smooth configurations play an important role. We may expect that a better result could 
be obtained once some other coefficient functions of the derivative expansion are added and the new set of coupled 
differential equations derived. 
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